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5— I ■ ABSTRACT 

Infrared Dark Clouds (IRDCs) arc cold, high mass surface density and high density structures, likely to 
be representative of the initial conditions for massive star and star cluster formation. CO emission from 
IRDCs has the potential to be useful for tracing their dynamics, but may be affected by depleted gas 
phase abundances due to freeze-out onto dust grains. Here we analyze C 18 J = 1 — > and J = 2 — > 1 
emission line data, taken with the IRAM 30m telescope, of the highly filamentary IRDC G035.39.- 
0^ , 0033. We derive the excitation temperature as a function of position and velocity, with typical values 

of ~ 7 K, and thus derive total mass surface densities, £ci80> assuming standard gas phase abundances 
and accounting for optical depth in the line, which can reach values of ~ 1. The mass surface densities 
reach values of ~ 0.07 gem -2 . We compare these results to the mass surface densities derived from 
^sO • mid-infrared (MIR) extinction mapping, Ssmf, by Butler & Tan, which are expected to be insensitive 

to the dust temperatures in the cloud. With a significance of > 10c, we find Eciso/^smf decreases 
by about a factor of 5 as S increases from ~ 0.02 to ~ 0.2 g cm' 2 , which we interpret as evidence for 
CO depletion. Several hundred solar masses are being affected, making this one of the most massive 
clouds in which CO depletion has been observed directly. We present a map of the depletion factor in 
the filament and discuss implications for the formation of the IRDC. 
Subject headings: ISM: clouds, dust, extinction — stars: formation 



1. INTRODUCTION 

Silhouetted against the Galactic background, Infrared Dark Clouds (IRDCs) are opaque at wavelengths ~ 10 /zm 
(Perault et al. 1996; Egan et al. 1998), cold (T < 20 K; Carey et al. 1998; Pillai et al. 2006), and dense (n H > 10 3 -10 5 cm 3 ; 
Teyssier et al. 2002; Rathborne et al. 2006; Butler & Tan 2009, hereafter BT09; Peretto & Fuller 2010). They are likely 
to be the precursors of massive stars and star clusters as they have similar physical conditions, such as mass surface 
densities, as regions with such star formation activity (Rathborne et al. 2006; Tan 2007; Zhang et al. 2009; Ragan et 
al. 2009). CO emission from these clouds may be useful for understanding their dynamics (e.g. Hernandez & Tan 2011, 
hereafter HT11), but could be affected by depleted gas phase abundances due to freeze-out onto dust grains, especially in 
the coldest, highest density regions. 

Gas phase depletion of CO, averaged along the line of sight, has been observed in the cold (T < 10 K) centers of 
relatively low-mass and nearby starless cores, (e.g. Willacy et al. 1998; Caselli et al. 1999; Kramer et al. 1999; Bergin 
et al. 2002; Whittet et al. 2010; Ford & Shirley 2011). Typically, depletion is characterized by measuring the depletion 
factor, fjj, defined as the ratio of CO column density expected assuming standard gas phase abundances given the column 
of material observed from either the mm dust continuum emission or near infrared (NIR) dust extinction to the observed 
CO column density (typically from C 17 or C 18 0). Caselli et al. (1999) estimated the expected CO column based on mm 
dust continuum emission, which has the advantage of being able to probe to high column densities, but is sensitive to the 
adopted dust temperature and emissivity. They concluded depletion affected a region at the core center containing about 
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2 Mq of gas, where Tin ^ 10 5 cm -3 , with depletion factors of up to ~10 where the mass surface density is £ ~ 0.6 g cm -2 . 
Kramer et al. (1999) estimated the expected CO column based on NIR extinction, which does not require knowing the 
dust temperature, but does require there to be a sufficient areal density of background stars detectable in the NIR. They 
found depletion factors of up to ~ 2.5 for regions with Ay ^20 — 30 mag, corresponding to £ ~ 0.1 — 0.15 g cm~ 2 . 

Massive protostellar cores and clumps are typically more distant and difficult to study, but CO depletion has been 
reported by Fontani et al. (2006) from a study of 10 sources with median fjj ~ 3.2 (but a dispersion of about a factor of 
10), Thomas & Fuller (2008) from a study of 10 sources with a mean fjj ~ 1.3 and Lo et al. (2011) from a study of 1 source 
with fjj ~ 10. These results rely on estimates of the expected CO column density based on mm dust continuum emission, 
are derived only for single pointings to the sources, and can depend on radiative transfer modeling of the unresolved 
source density and temperature structure (Thomas & Fuller 2008; Lo et al. 2011). Source to source comparisons are 
hampered by possible isotopic abundance variations affecting these rare CO isotopologues. The above sources already 
contain massive protostars, but it is not clear if the depletion signal arises from the immediate surrounding envelope or 
from nearby unresolved starless cores. Some of the massive protostars studied produce ultra-compact H II regions and 
photodissociation of molecules could be occurring in localized regions, which would mimic depletion. 

We expect CO depletion to be widespread in the dense regions of IRDCs, potentially affecting: the physical properties 
one derives from CO emission; the mid and far infrared opacities of dust grains as CO ice mantles build up; and thus 
the initial conditions of star and planet formation in these regions. Individual resolved IRDCs, assumed to have uniform 
isotopic abundances, may also be useful laboratories in which to study the depletion process as a function of local gas 
conditions. 

In this paper, we present IRAM 30m observations of C 18 J = 1 — > and J = 2 —> 1 emission from the filamentary 
IRDC G035. 30-00. 33 (Cloud H in BT09; near kinematic distance of d — 2.9 kpc). To look for evidence of depletion, the 
C 18 0-derived mass surface density, ScisOj is compared with the small median filter (SMF) mid-infrared (MIR) extinction 
mapping derived mass surface density, Esmf (BT09; Butler & Tan 2011, hereafter BT11). This work is motivated by 
the study of HT11, who used 13 CO molecular line emission from the Galactic Ring Survey (GRS) to estimate the mass 
surface densities of two highly filamentary IRDCs, including Filament H. Assuming a constant value of T cx = 15 K, HT11 
found tentative evidence for CO depletion, but could not exclude the possibility that other effects, such as systematic 
changes in the excitation temperature or the contribution of high opacity cores, were the cause of the observed decrease 
of Ei3Co/Esmf with increasing E. With our new higher-resolution, multi-transition C 18 data, we are able to exclude 
or mitigate these effects, as well as resolving higher mass surface density structures to probe a larger range of conditions 
where depletion may be occurring. 

2. MASS SURFACE DENSITY FROM MIR EXTINCTION MAPPING 

The 8 /jm SMF mass surface density, Esmf, map was derived at 2" resolution from the Spitzer IRAC band 4 (Galactic 
Legacy Mid-Plane Survey Extraordinaire [GLIMPSE]; Benjamin et al. 2003) image by comparing the observed intensity 
at each position with the expected background intensity, estimated by interpolating the intensities of surrounding nearby 
regions where median filter smoothing is used to define the background model (see Figure la and lb). Following BT09, a 
dust opacity of fi^m = 7.5 cm 2 g _1 was adopted, similar to the filter response and background spectrum weighted mean 
IRAC band 4 opacity expected from the Ossenkopf & Henning (1994) thin ice mantle moderately coagulated grain model 
with a gas-to-dust mass ratio of 156. This value is somewhat higher than values adopted by other dust models (e.g. 125 
is used for the Weingartner & Draine 2001), although a recent estimate from depletion studies finds a gas-to-dust ratio of 
141 (Draine 2011, p265). In any case, as described below, our study of CO depletion compares relative abundances as a 
function of £ in the IRDC and so is independent of this choice of overall normalization. 

A correction for foreground emission also needs to be estimated. BT09 made this correction by estimating the amount 
of foreground emission from a physical model of the Milky Way and given a measured kinematic distance (assumed to 
be near) of the cloud. Battersby et al. (2010) have pointed out an additional source of foreground from scattering in the 
IRAC array. BT11 have developed a more accurate empirical method for estimating the foreground emission, based on 
the presence of independent saturated (high optical depth) cores, and here we use this new method. For the region we 
analyze in this particular IRDC, the values of Esmf are increased by about 10% from those presented by BT09. 

Esmf in the filament is derived from comparison with adjacent regions, which are assumed to have negligible MIR 
extinction. In reality, we know from molecular line observations (e.g. 13 CO from the GRS analyzed by HT11), that 
these regions do have some material present associated with the IRDC. We refer to this as the IRDC "envelope" . The 
presence of the envelope and other systematic uncertainties associated with estimation of the MIR background intensity 
mean that Esmf becomes unreliable when < 0.01 gem -2 . For our comparison with the mass surface density derived 
from C 18 emission, the Esmf map is regridded to the much lower resolution of the CO data (see below) and all pixels 
with Esmf < 0.01 gem -2 are excluded from the analysis. Methods of accounting for the envelope material are discussed 
further in 

As noted by BT09, we must also account for locations of bright MIR emission. Wherever the observed MIR intensity 
is greater than the adopted background model an unphysical negative value of E will be estimated. Negative values of E 
are allowed up to levels comparable with the observed noise, but more extreme values, which are mostly due to discrete 
MIR bright sources, have Esmf set to zero. This causes an underestimation of the mass surface density in these regions. 
We identify and exclude from further analysis remaining (i.e. > 0.01 gem" 2 ) pixels in Esmf map (smoothed to the CO 
resolution) that have more than 20% of their area occupied by zero or negative values. In Figure lb, these excluded 
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pixels are indicated with "X" and "O" symbols for the CO(l-O) and CO(2-l) resolutions, respectively. Their exclusion is 
due either to the presence of a MIR bright source or in regions where the background modeling is inaccurate, which can 
sometimes occur near the edge of the filament. Only a relatively small number of pixels are affected by this exclusion. In 
fact, our final results would not have varied significantly if this exclusion had not been implemented. 

3. MASS SURFACE DENSITY FROM C 18 EMISSION 
3.1. Observations 

The C 18 J = 1 — > and J = 2 — > 1 lines were mapped using the IRAM (Institute de Radioastronomia Milimetrica) 
30m antenna in Pico Veleta, Spain in August and December 2008. An area of 2' x 4' was mapped using the On-The- 
Fly (OTF) method towards G035.39-00.33 with a central position of a(J2000) = 18 h 57 m 08 s , £(J2000) = 02°10'30" 
(/ = 35.517°, b = -0.274°). While the C 18 J = 1 -> transition was observed with the ABCD receivers with typical 
single side band (SSB) rejections > 10 dB, the C 18 J = 2 — >• 1 lines emission was mapped by using the HERA multi- 
beam receiver. Off-positions for both transition lines were set to (1830", 658"). The beam size at ~ 110 GHz for the 
J = 1 — > transitions is 22", while at ~ 220 GHz the J = 2 — > 1 beam size is 11" . The VESPA spectrometer provided 
spectral resolutions of 20kHz and 80kHz for the J = 1 — > and J = 2 — > 1 lines respectively, which correspond to 
velocity resolutions of ~ 0.05 km s _1 and ~ 0.1km s _ . For this study, all spectra were resampled to the same velocity 
resolution of 0.2 km s _ . The typical system temperatures were 150-220 K. Intensities were calibrated in units of antenna 
temperature (T^), and converted into a main beam brightness temperature, Tb,^, via Ta = ^/ciumpls.,,, where rj is 
a main beam efficiency and /clump is the beam dilution factor. We use rj = 0.64 for the J = 1 —> transition, and 
?7 = 0.52 for the J = 2 — >• 1 transition. The typical la RMS noise of the data is 0.2/^ km s _1 over the velocity range of 
40 — 50 km s . Since the C 18 emission is extended over the filament, we assume /dump = 1- Figure [TJ; presents the 
morphology of Filament H as seen in C 18 J = 2 — ► 1 emission. 

3.2. Mass Surface Density and T ox Estimates 

We estimate the column density of C 18 molecules, dAcisO) in the velocity interval dv, from their emission through 
the general equation: 

dN C iso(v) 8tt gi t v Q Iot 



dv AXl g u l- exp (-hv/kT CK ) giexp(-Ei/kT c ^) ' 



(1) 



Here Q ro t is the partition function for linear molecules given by Qrot = X)jLo(2^ + l)exp(— Ej / k T ox ) with Ej = J(J + 
l)hB, where J is the rotational quantum number and B is the C 18 rotational constant equal to 5.4891 x 10 10 s . 
hv/k = 5.269, 10.54 K for J = 1 — > and J = 2 — > 1 transitions, respectively. At 7.5 K, Q Iot — 3.205. A is the Einstein 
coefficient, 6.266,60.11 x 10 _8 s _1 for J = 1 — > and J = 2 — > 1, respectively. Ao is the wavelength of the transition, 
0.273, 0.137 cm for J = 1 — > and J = 2 — > 1, respectively, gi and g u are the statistical weights of the lower and upper 
levels, and t v is the optical depth of the line at frequency v, i.e. at velocity v. The excitation temperature, T cx , is assumed 
to be the same for all rotational levels. Details on the estimation of T ox are given below. 
The optical depth, r„, is derived through the detection equation: 

T B , v = ^[f(T ex )-f(T hs )][l-e-^] (2) 

where Tb d is the main beam brightness temperature at frequency v, f(T) = [exp(hv/[kT]) — 1] , and Xb g is the 
background temperature of 2.725 K. For the observable, Tb, V : and for an assumed T ex , t v can be solved for directly 
through equation |2]). Therefore, we can solve for the column density per unit velocity, dA'igco/di', at each I, b, v position. 

While care is taken to account for the optical depth in our column density estimates, for reference we also state the 
case of the optically thin limit of the C 18 (J = 1 —¥ 0) column density. If t v is small, then equation Q reduces to 
Tb,v = (hv '/ 'k)[j '(Tcx) — /(Tb g )]Tv. Inserting into equation (TTJ) gives: 

= 6.571 x 10" ^°) (T J l - exp(-W^c X )]-^ cm-km^s (3) 

dv /(T ex ) - /(Tbg) ?7/clump 

-)• 9.758 x 10 14 TaI*£ cm _2 km _1 s (T cx = 7.5 K). (4) 

T] Jclump 

As in HT11, an inspection of the C 18 emission in /, b, v space indicates that the gas associated with the filament is in 
the range of 40 — 50 km s _1 . The total column density per pixel is then calculated over the entire velocity range of the 
filament, N C i80 = jdNcmo- 

The column densities for both transitions, A^cisOj are converted to a total mass surface density EcisOj by assuming 
the abundance ratios of nieo/mso = 327 from Wilson & Rood (1994) and ni2Co/nH2 = 2 x 10~ 4 from Lacy et al. (1994). 
Thus, our assumed abundance ratio of C 18 to H 2 is 6.12 x 10~ 7 and E for each pixel is then given by: 

E clso = 7.652 x 10' 2 in ^ 18 ° 2 g cm- 2 , (5) 
10 ib cm~" ! 

assuming a mass per H nucleus of /in = 2.34 x 10~ 24 g, i.e. E — 1 g cm -2 is equivalent to A^h — 4.27 x 10 23 cnr 2 . 




Fig. 1. — Morphology and depletion maps of the IRDC. (a,) Top left: Spitzer GLIMPSE IRAC 8 (J>m image, with linear intensity scale in MJy Sr _1 . The image has 1.2" 
pixels and the PSF has a FWHM of 2". (b) Bottom left: Mass surface density, Esmf, with linear intensity scale in gem -2 , derived from the image in panel (a) using the 
small median filter (SMF) MIR extinction mapping method of Butler & Tan (2009; 2011). Regions with Esmf > 0.01 g cm -2 but which are > 20% affected by artifacts 
in the extinction map (e.g. due to MIR bright sources) are excluded from analysis and shown by "X"'s and "0"'s for CO(1-0) and (2-1) resolution grids, respectively, (c) 
Top middle left: Integrated intensity map of C 18 (J = 2 — y 1) emission over the velocity range of 40 — 50km s _1 , i.e. the gas associated with the IRDC (HT11), in linear 
units of K kms - and a pixel scale of 5". (d) Bottom middle left: The mean excitation temperature map weighted by the column density in K, with pixel size of 11". (e) 
Top middle right: Relative depletion factor (/£>) map for Case 1 (no CO envelope subtraction), (f) Bottom middle right: Relative depletion factor (f'r>) map for Case 1 
HiRes (no CO envelope subtraction, E derived at the CO(2-l) resolution), (g) Top right: Relative depletion factor map for Case 2 (CO envelope contribution estimated 
via interpolation across strips 2, 3 and 4 then subtracted; note we consider this process unreliable for strip 1). (h) Bottom right: Relative depletion factor map for Case 2 
HiRes. 
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In order to accurately derive the mass surface density of the filament, an estimate of the excitation temperature, T ex , 
is needed. To perform this estimate throughout the filament, we varied the assumed temperature at each l,b,v position 
until the ratio between the column densities derived from both transitions were in agreement. To do this, we first defined 
i?2,i as the ratio between the J = 2 — > 1 and J = 1 — > column densities: 

_ d7V C 180,21 

#2,1 = TT7 • (6) 

aJv C i 80 ,io 

This method is similar to the one used in Kramer et al. (1999), except they averaged over the velocity profile of their 
cloud. The higher resolution J = 2 — » 1 data was convolved with a beam of 22" and regridded to match the resolution and 
pixel scale of the J = 1 — > data. For all /, b, v positions above a noise limit of 3er in both transitions, i?2,i was calculated 
first assuming a T cx = 30 K. Then, T ex was iteratively decreased until i?2,i converged to unity. This step provided a three 
dimensional grid containing estimates of T ex for all positions above the noise limit. Next, for all positions below the noise 
threshold, their T ex was estimated by taking the mean excitation temperature at the corresponding I, b position. Finally, 
for any remaining l,b, v positions without an estimated excitation temperature, the mean T ex of 7.2 K resulting from the 
previous steps was used. Positions left for this final step are mainly in the outer regions of the filament where the emission 
is weak and/or the noise is high. The column density weighted T ox map is shown in Figure [T}1. 

4. COMPARISON OF S C 180 AND Ssmf: 
EVIDENCE FOR CO DEPLETION 

In Figure [TJ we present the morphology of the filamentary IRDC H. The goal of this section is to compare Sciso and 
Ssmf- The simplest way of doing this, which we refer to as Case 1, involves a straightforward pixel by pixel comparison 
of these values, smoothing the Ssmf data to the resolution of the C 18 O(T0) observations, for which we have derived 
accurate excitation temperature information. Note, that only pixels with Ssmf and Sciso > 0.01 g cm~ 2 are considered. 
Also, pixels for which Ssmf is affected by bright MIR emission are excluded (see $Z\). We also perform a comparison 
at the higher angular resolution of the C 18 0(2-l) observations, which we refer to as Case 1 HiRcs, assuming T cx at this 
higher angular resolution can be estimated from the values derived at the (1-0) resolution. For both these versions of 
Case 1, we refer to Sciso as Sciso.tot, since it is derived from all the C 18 emission associated with the IRDC and its 
surrounding GMC. 

However, as is apparent from Figure [TJ the C 18 emission is more extended than the 8^m extinction map from Butler 
& Tan (2009; 2011). This is because, as discussed above, the extinction map is derived from an "on-off" comparison 
with adjacent regions, which help define the background MIR intensity that is expected to be behind the filament. Thus 
the MIR extinction mapping method becomes insensitive to material present in these adjacent, lower column density 
("envelope") regions. A fair comparison between Ssmf and Sciso would allow for this envelope material. We thus define 
"filament" and "envelope" regions based on the 8 /im image of the IRDC. Following HT11, the filament is defined to be 
a rectangular strip centered at a(J2000) = 18 h 57 m 08.02 s , <5(J2000) = 02°10'35.7", 2.05' wide in R.A. and 4.47' long in 
Dec. The outline of this filament region is shown by a red box in the panels of Figure 1. The envelope region is defined 
to be made up of two adjacent rectangular regions on either side of the filament. These are shown as blue rectangles in 
Figure 1 and are each 0.56' wide in R.A. and 4.47' long in Dec. Note, that because of the limited area mapped by our 
observations, these envelope regions are narrower than those considered by HT11. 

For our Case 2, we assume that the C 18 material present in the envelope regions is also present at the similar levels 
towards the filament region, and so attempt to subtract this emission from the C 18 spectrum of the filament, before then 
comparing to Ssmf- To carry out this subtraction we divide the filament and envelope into four E-W strips (1 to 4 from 
N to S) (see Figure 1). In each strip, the mean column density per unit velocity is evaluated for the filament (based on 
66 C 18 O(l-0) pixels) and the two adjacent envelope regions (based on 18 C 18 O(l-0) pixels each) (see Figure[2]), using the 
Tex estimates described previously. The envelope spectra are averaged and then subtracted from the filament. The total 
column of this envelope-subtracted spectrum is evaluated and used to derive Sciso,fil- This is of course an approximate 
method for accounting for the envelope material: one can see from Figure [2] that the envelope spectra on either side of 
the filament can be quite different, especially for strips 1 and 2. The uncertainty in the envelope-subtracted spectrum 
becomes large when the envelope spectra are of similar strength as that of the filament, as is the case for strip 1. Thus we 
do not regard the results of envelope subtraction for strip 1 as being reliable, and we exclude these pixels from the Case 
2 analysis. As with Case 1, we also perform a Case 2 HiRes analysis, using Sciso,fil estimated at the higher resolution 
of the CO(2-l) data, adopting values of T cx evaluated at the CO(1-0) resolution. 

With these Case 1 and 2 methods, we now compare the pixel by pixel values of Eciso with Ssmf derived from 
MIR extinction mapping. As noted in HT11, these measurements of S are essentially independent of cloud distance 
uncertainties. Figure presents Sciso,tot versus Ssmf, i-e. Case 1 of no envelope subtraction. The best fit power 
law relation to the CO(1-0) resolution data of Sciso, Tor/g cm~ 2 = A(£sMF/g cm -2 )" has a = 0.452 ± 0.054 and 
A = 0.146 ± 0.023. For Case 1 HiRes (i.e. at the CO(2-l) resolution, adopting CO(1-0) resolution T GX estimates) we find 
a = 0.463 ± 0.025 and A = 0.151 ± 0.010. These results are summarized in Table [JJ 

These uncertainties are derived assuming that the errors of each individual measurement are as follows: for Sciso, a 
fixed value of 0.0024 gem -2 (derived from the ler RMS noise of 0.2Kkm s _1 over the velocity range of 40 — 50km s _1 ) 
and a 20% error to account for uncertainties in T cx assumed to be 1 K at the typical temperature of 7 K; for Ssmf, a 
15% error plus a systematic error of O.Olgcm -2 (BT09). At the resolution of the CO pixels (ll"for CO(1-0) and 5"for 
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C0(2-l)), the Esmf measurements are independent, but the Eciso results are not since the telescope beam is about twice 
the pixel scale. Thus the above quoted uncertainties of the power law fits assume, conservatively, only 25% of the pixels 
are used (although the derived values of the parameters are based on fits to all of the pixels). 

We argue below that Esmf is a more accurate measure of the true mass surface density in IRDCs than Eciso j since 
one does not expect large changes in MIR dust opacities in these environments, based on the Ossenkopf & Hcnning (1994) 
dust models. If this is true, then if C 18 were also an accurate tracer of mass surface density, then we should see a 
one-to-one relation between Eciso and Esmf, i-e. a~l, even if A (the value of Sciso/^SMF when Esmf — 1 gem -2 ) 
is not exactly unity because of systematic uncertainties in the absolute values of C 18 abundance or MIR dust opacities. 
We measure a = 0.452 ± 0.054 for Case 1 and a = 0.463 ± 0.025 for Case 1 HiRes, which are significantly (IOct and 21<r) 
different from one, and we interpret these results as being evidence for CO depletion from the gas phase. 

To illustrate that these results do not depend on the choice of dust opacity per unit gas mass, we have repeated the 
analysis but with a gas-to-dust mass ratio of 100 (rather than our fiducial value of 156). We find a = 0.509 ±0.073 (about 
7<T different from a = 1) for Case 1 and a = 0.552 ± 0.035 (about 13cr different from a = 1) for Case 1 HiRes. Note that 
we do not expect to derive exactly the same values of a as before since we have a fixed threshold of E > 0.01 gem -2 
to include points in the analysis and so reducing the gas-to-dust mass ratio causes us to lose some data points near this 
limit. 

Figure shows the ratio Eciso, tot/Esmf versus Esmf for our fiducial Case 1 and Case 1 HiRes analyses, with the 
derived power law relations overlaid. For 0.01 < EsMF/gcm~ 2 < 0.03 the mean values of Eciso, tot/Esmf are 1-316 
and 1.471 for Case 1 and Case 1 HiRes, respectively. By the time Esmf ^ 0.1 gem -2 , Eciso, tot/Esmf has declined to 
values of < 0.4. 

In Case 2 we attempt to account for the IRDC envelope: we consider that we can do this reliably only for strips 2, 
3 and 4, where the envelope is relatively weak compared to the filament. Figure presents Eciso, fil versus Esmf for 
Case 2. The best fit power law relation to the CO(1-0) resolution data of Eciso, FiL/g cm~ 2 = yl(EsMF/g cm -2 )" has 
a = 0.239±0.080 and A = 0.074±0.017. For Case 2 HiRes (i.e. at the CO(2-l) resolution, adopting CO(1-0) T ex estimates) 
we find a = 0.317 ±0.038 and A = 0.090 ±0.010. These uncertainties assume the same measurement uncertainties as Case 
1, except an additional systematic error of 0.01gcm~ 2 has been applied to Eciso, fil °-ue to uncertainties associated with 
envelope subtraction. Again these results indicate a significant (lOcr and 18cr for Case 2 and Case 2 HiRes, respectively) 
departure from a one-to-one (a = 1) relation, which we again interpret as evidence for CO depletion. The results with 
a gas-to-dust mass ratio of 100 are a = 0.303 ± 0.11 (about 6a different from a = 1) for Case 2 and a — 0.372 ± 0.048 
(about 13(7 different from a = 1) for Case 2 HiRes. 

Figure[3ji shows the ratio Eciso fil/Esmf versus Esmf for Case 2 and Case 2 HiRes, with the above power law relations 
overlaid. For 0.01 < EsMF/gcm~ 2 < 0.03 the mean values of Sciso, fil/Esmf are 1-099 and 1.238 for Case 2 and Case 
2 HiRes, respectively. These values are smaller than their equivalents for Case 1, as is to be expected now that we are 
allowing for the molecular envelope. The values are also very close to unity, suggesting that our adopted C 18 abundances 
and dust opacity per unit gas mass are reasonable. Again, by the time Esmf ^ 0.1 g cm -2 , Eciso, tot/Esmf has declined 
to values of < 0.4. 

4.1. Alternatives to CO Depletion 

There are several physical processes that could be responsible for the observed trend of decreasing Eciso /Esmf with 
increasing Esmf- One possibility could be that our corrections for the optical depth of the C 18 emission are systematically 
underestimated near the center of the filament where the column density is large. However, the largest optical depth 
corrections in the highest column density locations increase the column by only 30% (the highest optical depths are ~ 1, 
but lower when averaged over the whole column), so this effect is unlikely to be driving the observed trend. 

HT11 suggested their observed trend of decreasing Ei 3 co/^SMF with increasing Esmf could potentially result if at the 
same time there is a systematic decrease in the excitation temperature of about 5 K. However, from our T ex estimates, 
we find no strong negative temperature gradient within the IRDC towards the mass surface density peaks. In fact, T cx 
increases slightly towards to the center of the filament, probably as the densities become greater than the effective critical 
densities and the lower CO levels can thermalize. Thus, we exclude trends in T ex as causing the observed variation of 

Sci8o/^SMF- 

Fractionation of C 18 could in principle change the local abundance of this molecule, but the most important way in 
which this can be achieved is via isotope selective photodissociation at cloud edges, which would not be able to explain 
the trends of decreasing C 18 abundance that we see running from E ~ 0.02 gem -2 (Ay — 4 mag) to ~ 0.2 gcm~ 2 
(A v ~ 40 mag). 

Another possibility to be considered is systematic changes in 8 /im dust opacities for gas at higher densities. If the 
opacity was to increase (e.g. due to grain coagulation and/or ice mantle formation and growth), then this could explain 
our observed trend of decreasing Eciso/Esmf with increasing E. The Ossenkopf & Henning (1994) dust models do show 
an increase of Ks^m of 19% going from the uncoagulated thin ice mantle model to the uncoagulated thick ice mantle (all 
volatiles depleted) model. Maximal coagulation (corresponding to that expected after 10 5 yr at densities of 10 s cm~ 3 or 
after ~ 10 8 yr at densities of ~ 10 5 cm -3 , which is probably more that can be expected to have occurred since the 
observed densities of IRDC cores are < 10 5 cm -3 ; BT09) raises res/im by an additional 17%. Thus, ice mantle growth and 
grain coagulation appears to be able to account for only a small fraction of the observed variation of Eciso/^smf- 

We conclude the most likely cause of the trend of decreasing Eciso /Esmf with increasing Esmf is CO depletion due 
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to freeze out onto dust grains. This would cause a systematic reduction in the amount of CO gas observed in higher mass 
surface density regions, which are likely to also be of higher volume density. 



where X^ Q is the expected abundance of CO relative to H 2 given standard gas phase abundances, X° Q is the observed 
abundance and the last equality assumes that Esmf estimated from MIR extinction mapping is an accurate measure of 
the true mass surface density (this assumption is discussed further below). Given the uncertainties in the absolute values 
of the C 18 abundance and the MIR dust opacity per unit gas mass, we renormalize fr> to be unity for the regions of the 
IRDC with 0.01 < £sMF/gcm~ 2 < 0.03 and refer to this renormalized value as the relative depletion factor f' D = Bfp, 
where the scaling factor, B = 1.316, 1.471, 1.099, 1.238 for Case 1, Case 1 HiRes, Case 2, Case 2 HiRes, respectively. We 
show maps of f' D for these four cases in Figure [T^-h. We note that the values of f' D presented here, peaking at values 
~ 5, are mass surface density weighted averages and thus lower limits to the maximum values of the depletion factor that 
apply in the densest regions of the cloud. 

We conclude that with high (~ lOc) significance, widespread CO depletion is occurring in this IRDC, with depletion 
factors of up to ~ 5 (see Table [T]). These values are larger than those seen towards more evolved cores and clumps already 
containing massive protostars (Fontani et al. 2006; Thomas & Fuller 2008). Our measurement of CO depletion suffers 
from fewer systematic uncertainties, especially since we do not require knowledge of the dust temperature. 

Each pixel in the lower resolution depletion maps (11", half the C 18 O(l-0) angular resolution) corresponds to a length 
of 0.155 pc at the cloud distance of 2.9 kpc, and so contains a mass of 11.4(E/0.1g cm~ 2 ) M©. Thus, hundreds of solar 
masses appear to be affected by depletion along the filament (the total SMF-derived mass in the 4 strips is 580 ± 230 M , 
HT11), including a particularly prominent massive core or clump in strip 2 and a larger clump partially in strip 4 and 
extending to the south. 

Thus, IRDC G035. 30-00. 33 is one of the most massive clouds in which CO depletion has been detected by direct CO- 
based and non-CO-based measurements of mass surface density. Our results also suggest that CO depletion will be a 
common occurrence in IRDCs, since the values of E ~ 0.1 gem -2 in this cloud are quite typical (e.g. BT09). CO is 
therefore an imperfect tracer of a significant fraction of the mass of IRDCs (not just the coldest, densest cores). Accurate 
accounting for depletion and/or use of species suffering minimal depletion, such as NH3 and N2H + , are required for more 
accurate dynamical studies of these clouds. 

An estimate of the CO depletion timescale due to freeze-out onto dust grains is to — 8000/(nH 2 ,5<5) yr, where nn 2 ,5 is 
the number density of H2 molecules in units of 10 5 cm -3 and S is the sticking probability (of order unity; e.g. Tielens & 
Allamandola 1987) for CO on grains. We can apply this to the thinnest region of the IRDC: the ~ 5" (0.070 pc) wide 
filament near the center of strip 3, which appears to have significant CO depletion with f' D ~ 3 — 4. Assuming the depth 
of the filament, which has Esmf — 0.2 gem -2 , is similar to its width, then riH 2 ,5 = 2.0 and ~ 4000 yr. This provides 
a lower limit to the age of this part of the IRDC. The free-fall time, tg = (37r/[32Gp]) 1 / 2 , for this density is 6.9 x 10 4 yr, 
i.e. much longer. However, if the filament has been created by larger scale supersonic flows, then one might expect the 
high density gas to have been present for about the flow crossing time across the width of the filament. Velocities of 
~ 10 km s _1 may be relevant in models of GMC-GMC collisions (Tan 2000) or if the large-scale SiO emission seen towards 
this filament (Jimenez-Serra et al. 2010) has been created by such flows. The flow crossing time at this speed for this 
part of the IRDC is only 6800 yr. Thus the fact that we see CO depletion in these very thin filaments of the IRDC can 
help to constrain models for the cloud's formation. For models in which the cloud lifetime is less than the flow crossing 
time across the filament, a constraint is placed on the flow speed. For the thinnest region of this IRDC, this corresponds 
to flow speeds < 17 km s _1 . 

We thank E. van Dishoeck and R. Visser for helpful discussions and the comments of an anonymous referee, which 
helped improve the paper. AKH acknowledges support from a SEAGEP Dissertation Fellowship. JCT acknowledges 
support from NSF CAREER grant AST-0645412; NASA Astrophysics Theory and Fundamental Physics grant ATP09- 
0094; NASA Astrophysics Data Analysis Program ADAP10-0110 and a Faculty Enhancement Opportunity grant from 
the University of Florida. 
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Table 1 



Parameters of Depletion Factor Analysis 



Case 



a 



A 



B 



/p(max) 



Case 1 

Case 1 HiRes 
Case 2 

Case 2 HiRes 



0.452 ±0.054 
0.463 ±0.025 
0.239 ±0.080 
0.317 ±0.038 



0.146 ±0.023 
0.151 ±0.010 
0.074 ±0.017 
0.090 ±0.010 



1.316 
1.471 
1.099 
1.238 



3.5 (at E S mf 

4.6 (at Esmf 

3.8 (at Esmf 

4.9 (at Esmf 



0.16 gem 2 ) 
0.20 gem" 2 ) 
0.16 g cm -2 ) 
0.20 gem" 2 ) 
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Fig. 2. — Velocity structure of the C O molecules associated with the IRDC and its envelope. The column density distribution, 
diVc^o/di 1 , has been derived from the C 18 O(l-0) and (2-1) spectra, local estimates of T cx and including optical depth corrections. 
(a.) Left: the 4 sets of profiles (offset to display from top to bottom and labeled 1 to 4) correspond to the 4 strips shown in 
Figure [T] The dotted, red line is the summed contribution from gas from the central region of each strip, corresponding to the 
IRDC "filament" (see Figure Q] and text). The dot-dashed and long-dashed blue lines show summed contribution from the gas 
from the eastern and western envelope regions, respectively, (b) Right: Illustration of envelope subtraction (Case 2, see text). For 
the same strips as in (a), we subtract the average of the eastern and western envelopes (short-dashed blue lines) from the filament 
(dotted red lines), to leave an estimate of the material in the filament (solid black lines). We consider this process unreliable for 
strip 1, where the envelope contains a similar amount of material as the filament. 
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Fig. 3. — Evidence for CO depletion, (a) Top Left: Comparison of E c iso,tot (i- e - Case 1) and Esmf for all C 18 O(l-0) (crosses) 
and C ls O(2-l) (dots) pixels for which both Eciso and Esmf > 0.01 g cm -2 and the pixel is < 20% affected by Esmf artifacts, e.g. 
due to MIR bright sources. The dotted line shows the condition E c iso tot = Esmf- The solid, dashed lines show the best-fit power 
law relations to the C 18 O(l-0), C ls O(2-l) resolution data, respectively, (b) Bottom Left: Ratio Eciso tot/Esmf (i.e. Case 1) 
versus Esmf, with the same symbol and line notation as in (a). The horizontal solid, dashed lines from 0.01 < Esmf /g cm" 2 < 0.03 
indicate the mean values of the data in this range for the C 18 O(l-0), C 18 0(2-l) resolution data, respectively. The cross in the 
upper-right corner indicates typical estimated uncertainties, (c) Top Right: Same as (a), but now estimating E c i8 FIL from 
molecular gas associated with the filament after envelope subtraction (Case 2) in strips 2, 3 and 4. (d) Bottom Right: Same as (b), 
but for Case 2. Both (b) and (d) show that Eciso /Esmf decreases by up to a factor of ~ 5 as Esmf increases from ~ 0.02 g cm" 2 
up to ~ 0.2 g cm -2 . 



